Lateral Effects in Fermion Antibunching 
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Lateral effects are analyzed in the antibunching of a beam of free non-interacting fermions. The 
emission of particles from a source is dynamically described in a 3D full quantum field-theoretical 
framework. The size of the source and the detectors, as well as the temperature of the source are 
taken into account and the behavior of the visibility is scrutinized as a function of these parameters. 
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I. INTRODUCTION 



The wave function of two fermions is antisymmetric under exchange of the two particles, as a consequence of the 
Pauli exclusion principle. For this reason, the probability amplitude for their being spatially close together is small 
and their correlated detections are reduced when compared to a random sequence of classical particles. This very 
distinctive quantum feature is named antibunching and has no classical analog. Notice that in general the two particles 
can be emitted from totally incoherent sources. 

The analogous phenomenon for bosons is a cornerstone in the study of quantum correlations and was hrst observed 
in astronomy, where it is known as the Hanbury Brown-Twiss effect Photon second-order coherence effects 
@j Hi EIj yielding bunching, are discussed in physics textbooks 0, 0], and led to novel interesting applications in 
quantum imaging 0, Q and lithography . 

The most relevant difference between the Bose-Einstein and Fermi-Dirac statistics are the phase space densities 
(occupation numbers), that change by several orders of magnitude. In a laser beam, one obtains values of order 10 14 , 
, while typical densities for thermal light, synchrotron radiation and electrons are of order 10~ 3 ; finally, for the most 
advanced neutron sources, one gets 10 -15 . These figures make it very difficult to observe fermion antibunching. In 
, addition, for charged particles (electrons and pions), additional Coulomb repulsion effects should be considered, that 
tend to reduce the visibility and mask the observation of the phenomenon. 



Quantum correlations have been detected in a series of interesting experiments: in condensed-matt er p hysics, 



where the electronic states are confined within the Fermi surface |h| El fl2| . for superconductor emitters [l3j . in the 



coincidence spectrum of neutrons from compound-nuclear reactions at small relative momentum [lil [l5j , as well as 
in pion pairs emitted from a quark-gluon plasma fl6| . Recently, antibunching was observed on a beam of thermal 
neutrons emitted from a nuclear reactor [ijj . This can be considered as a direct experimental evidence of free fermion 
antibunching, in which an ensemble of free Fermi particles displays quantum coherence effects. Other remarkable 
1 antibunching experiments have been recently reported for neutral atoms, both in a degenerate atomic Fermi gas (l8j 
5^ ! and in Fermi/Bose gases [l9fl . 

Huge numerical differences in phase space densities, like the afore-mentioned ones, call for close scrutiny of the 
theoretical premises as well as dedicated experimental efforts. Notice that these quantum statistical effects appear- 
to play a prominent role in phenomena that are characterized by figures that differ by almost 30 (!) orders of 
magnitude. The present study is motivated by this observation. We intend to analyze the antibunching phenomenon 
in the correlated detections of two neutral fermions, such as neutrons, emitted by a generic thermal source at a given 
temperature. Notice that bunching effects from (pseudo-)thermal sources still raise controversial interpretations (2(i| 
and are therefore worth investigating from first principles. 

Our main objective will be to analyze the spatial coherence and in particular the coherence area and volume. Lateral 
effects are becoming a critical issue, in view of a new generation of experiments. They were carefully analyzed in a 
series of experimental articles on X-ray bunching [2j|. For the sake of concreteness, we will focus on fermions, but 
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our analysis can be very easily extended to bosons (by replacing the energy distributions of the source and changing 
relevant signs in the formulas). 

We will treat both the source (a thermal oven) and the particle beam as fully (second) quantized systems and 
will study the emission process at thermal equilibrium, when the beam has reached its stationary configuration. 
This approach will have the advantage of treating both the oven and the fermion beam on an equal footing and of 
introducing the properties of the source in a natural way. 



Before starting a detailed analysis, let us outline the main features of the setup we have in mind and stress the 
main points of our argument. Our setup is the simple one schematically shown in Fig. [TJ Particles are emitted from 
a source through a small window, go through a monochromator (not shown), and are detected by two detectors. We 
count the number of coincident detections. At the initial time t = 0, the source is in the thermal equilibrium state at 
a finite temperature and outside there is the vacuum. Starting from this initial condition, we shall solve the dynamics 
of the emission, so that a stationary beam of particles will be prepared at t — > oo, after a transient period. The beam 
profile will not be added "by hand," but will be obtained by solving the equations of motion, so that the coherence 
properties of the emitted particles will reflect the dynamics of the emission. 

The lateral features of the system affect the antibunching, even when both detectors are placed on the longitudinal 
axis and we shall look at the correlation in the longitudinal direction. To this end, the lateral size of the detector 
mouth must be duly taken into account. We shall therefore implement the lateral resolution of the detector, as well 
as the longitudinal one, in the two-particle distribution function. The variables and parameters that characterize the 
setup are summarized in Table [TJ 

We shall start by writing down the Hamiltonian of this many-body system in Sec. IIII1 This is the crucial part of 
the present analysis, since it fully relies on dynamical consideration. In order to facilitate the introduction of the 
characteristics of the source, like temperature, size of the window, and so on, in a natural way, we shall adopt a 
two-field approach: one field describes the particles in the source and the other one the emitted particles outside. The 
emission Hamiltonian Aif cm i ss i on (which converts a particle in the source into a particle outside and vice versa) is at 
the heart of our analysis and must fully take into account all important features of the experimental setup, as well 
as the main characteristics of the physics of the emission process. The Hamiltonian below will enable us to discuss 
the lateral coherence features of the emitted beam, yet it will be simple enough to be (almost) solvable. As we will 
see, the diffraction of the particles emitted through the window governs the lateral coherence and is controlled by the 
lateral size of the emitting window. Once the Hamiltonian is written down, one has "only" to solve the equations of 
motion (and has no "freedom" anymore). 

The article is organized as follows. The dynamics of the emission is perturbatively solved in Sec. IIV1 under the 
assumption of weak emissivity, namely weak coupling A < 1, and the stationary limit t — > oo realizes a nonequilibrium 
steady state. The beam profile thus prepared is studied at a large distance from the source in Sec.[V] We then compute 
the two-particle distribution function, or in other words, the second-order correlation function, defined in Sec. I VII The 
interplay between the singlet and triplet contributions determines to which extent the coincidence counts are reduced 
(antibunching) when the two detectors are close to each other. Indeed, the singlet contribution yields bunching and 
the triplet one antibunching, with the latter three times larger than the former. The detector sizes (resolutions) a and 
d are implemented into the correlation functions, the saddle-point approximation is carefully worked out for the case 
in which the detectors are placed on the longitudinal z axis, and we obtain a formula for the normalized two-particle 
distribution function. The noncollinear case, with the two detectors placed off the longitudinal axis, is also discussed. 

The antibunching is then discussed and the coherence properties are clarified in Sec. lVIIl on the basis of our formula 



II. SETUP AND OUTLINE 




FIG. 1: Coincidence between two detectors in the beam of emitted particles: the interference of the two alternatives yields 
antibunching. 
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TABLE I: Summary of the variables and parameters used in the calculation. 



z longitudinal direction 

x, y transverse directions 

w lateral size of the circular emitting window of the source 

w z depth of the emitting region 

feo average momentum at monochromator 

Ski (i =-L,z) monochromator window 

a lateral size of the circular mouth of the detector 

d detector resolution in the longitudinal direction 

(origin) center of the emitting region 

fi (i = 1, 2) centers of detector apertures 

(3 inverse temperature of the source 

/i Fermi level (in the source) 

g(r) emitting window function 

f(k) monochromator momentum-window function 

Rr(r) detector resolution function 



for the collincar case, and the effect of the temperature of the source is scrutinized. The temperature effect is shown 
to be very weak. The dependence of the antibunching correlation function on the distance between the two detectors 
is found to be controlled by the lateral monochromator window and the longitudinal detector resolution, while the 
magnitude of the antibunching effect is determined by the lateral size of the source. Finally, a variety of experiments 
are analyzed in Sec. IVIII| in the light of the lateral coherence, and the main results are summarized in Sec. IfXI 



III. HAMILTONIAN AND STATE OF THE SOURCE 



Let us start with the Hamiltonian: we take 

H = Hq + Ai?emission 7 -^0 = ^bcam + ^source, (3.1a) 




Not only the emitted particles but also the source is treated as a second quantized dynamical system. The Hamiltonian 
of the particles in the source is H SOUTCe and that of the emitted particles is -ffbeam- The emission is dynamically described 
by an emission Hamiltonian A-ff e mission : a particle of momentum hi and spin a is annihilated by aw a in the source 
and is created by cj, CT outside with an amplitude XT^k' , the meaning of which will be described below. The creation 
and annihilation operators obey the canonical anticommutation relations for fermions 

i c k*,c{>a>} = W^lvi = 5 aa ,5 3 (k - k'), {c feCT ,a^ CT ,} = 0. (3.2) 

It is assumed that no spin flip occurs during the emission process and that the emission is irrespective of the spin 
state of the particle (TW does not depend on a): generalizations to more general cases are straightforward. The field 
operator in the configuration space for the emitted particles is denoted by 

Mr ' t)= J^ Ck ° el(k ' r ~ eht) - (3 - 3) 
In the following discussion, the dispersion relations are assumed to be = uik = k 2 /2m. 



4 



In Eq. (|3.ip . A is a small parameter, that will enable us to work in the weak-coupling limit. Although this approach 
is familiar in variety of theoretical approaches aimed at explaining diverse experimental situations, a few words of 
explanation are necessary in this case. We have in mind a situation in which an oven emits a beam of particles through 
a small aperture (which we refer to as "source"). Usually those particles that leave the source are monochromatized 
and can travel in waveguides, undergoing all kinds of losses. The parameter A globally accounts for all these diverse 
processes and A-ff mission simply enables us to take a particle with approximately the right characteristics in the 
oven and put it in the final section of the beam. The smallncss of the opening and the total "efficiency" of the 
emission process (from the oven to the region of space where the experiment is practically done, passing through 
monochromators, optical elements and/or waveguides and undergoing losses) calls for an approach in which A is a 
small parameter. We anticipate that in all final formulas, where normalized distribution functions will be studied, A will 
always simplify, making the final results independent of the details of the apparatus (such as the monochromatization 
procedure, reflection and transmission processes, losses in optical elements and waveguides and so on). Of course, one 
must be able to retain all essential elements in the analysis and final formulas. The quantity T kk ' in Eq. (|3.1c[) takes 
into account the action of the monochromator and the size of the source, and will be defined in Eq. (|3.5p . 

The Hamiltonian discussed in this section is to be considered as a phcnomenological transfer Hamiltonian, con- 
veniently tailored in order to discuss lateral size effects. It is similar to a "tunneling" Hamiltonian (for a two-field 
formulation of a tunneling process, see (22J) and can describe the particle emission from a small opening. 



A. Emission 



We consider the following emission process. Only the particles around the window of the source are emitted outside. 
That is, a particle in the momentum state |fc) (with k z > 0) is annihilated by a ka around the window of the source 
and is converted into a particle outside by c k , a . The emitting region is specified by a function g(r) centered around 
the window of the source that characterizes the lateral size of the window. One may further put a monochromator 
f(k) after the emission. The emission Hamiltonian is then given by 

Emission = £ f d 3 k' f d 3 k c{, a (k' \ f(k)g(r )6{k z ) I k) a ka + h.c. (3.4) 

with the emission matrix 

T k , k = (k'\f(k)g(r)e(k z )\k) = f{k')B{k z ) J ^ g(r)e-^ k ' = f(k')g(k' k)9{k z ). (3.5) 

The theta function, 8{k z ) = 1 for k z > and 0(k z ) = for k z < 0, accounts for the positivity of the longitudinal mo- 
mentum k z . In the following calculation, we assume Gaussian shapes for the emitting region and the monochromator, 



g(r)= . e -r-w-y2 W 2 

^(27r) 3 dctW 2 



f f k ) = t _ p -(k~ k0 ). {S K)-i( k - k0 )/± , 6]C y 

^(2^) 3 dct(<5/C) 2 



(w 2 

w 2 | 
^ i 

( (Sk^) 2 
(Sk ± ) 2 

V o o (Sk z ) 2 



where w represents the lateral size of the window of the source, w z is the depth of the emitting region, and Ski 
(i =_L, z) characterize the monochromator. In the following, we shall take fco = (0, 0, ko). 

It is interesting to notice that the non-factorized form (|3.5| of the interaction Hamiltonian will produce the required 
diffraction effect. (A factorized emission Hamiltonian would correspond to a point source, irrespectively of the state 
before the emission, and would not yield the desired lateral effect. The choice of the Hamiltonian and the validity of 
our working assumptions will be continuously checked throughout the whole calculation.) A particle with momentum 
k in the source is converted into a particle with momentum k' outside. The momentum transfer is governed by the 
Fourier transform g(k' — k) of the "interface" function g(r) and is ruled by the size of the window. A smaller window 
yields larger momentum transfer and results in a larger divergence of the emitted beam. This point will be crucial 
in the following discussion on the lateral coherence of the beam. The beam profile is therefore a direct consequence 
of the dynamics and is not artificially imposed at the outset. We also notice that the longitudinal component of 
momentum is not necessarily preserved during the emission process, as conservation of the longitudinal momentum 
prevents beam divergence. This motivates the choice of the form factor. 
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B. State of the Source 

Having written the Hamiltonian of the emission process, it is straightforward to introduce also the properties of the 
source. The initial thermal state of the source at a finite temperature /3 _1 is characterized by 

(«iUfcv<) = N(u; k )S aa ,5 3 (k - fc'), N(w) = - 1 . (3.8) 

In the present article, we shall focus for concreteness on the Fermi distribution. However, we can think of a more 
general distribution N(u>k). In fact, many of the formulas below remain valid as long as the initial state is stationary 
with respect to Ho, admits the Wick decomposition, and N(ujk) is a slowly varying function around ko- 

IV. DYNAMICS OF EMISSION 

The Hcisenberg equations of motion read 

i-^Cfe CT (t) = e k c kcr {t) + X f d 3 k' T kk *a k < a (t), 

d r . (41) 

^a ka {t) = u k a ka (i) + A / d 3 k' T k , k c k > a (t). 
By formally integrating the second equation and inserting it into the first, we get the equation for c ka (t), 

i^c k At) = e k c kr7 (t) + X J d 3 k / T kk , C - luJ «' t a k , r7 -i\ 2 dt' j d 3 fe' K kk ,{t - t')c k , a (t'), (4.2) 

where 

K kk ,{t)= [d s k"T kk „e- lu >>" t TZ, k „. (4.3) 
The intcgro-differential equation (|4.2[) is conveniently solved by Laplace transformation and the solution is given by 

c k „{t) = J d 3 k'G kk ,(t)c k , a -i\J^ dt' J d 3 k' J d 3 k" G kk .(t - t'jTv^e-"-*"*' a k „ a , (4.4) 

where 

GW(t) = [ ^G kk ,(s)c st , (4.5a) 

G k U s ) = ( s + 1£ k)S 3 (k - k') + A 2 /^(,s), K kk ,(s) = f d 3 k" Tkk " T *'k" , (4.5b) 

J S + lWfc" 

with Cb running parallel to the s-imaginary axis (Bromwich path). 

A. Nonequilibrium Steady State 

To take the stationary limit t — ^ oo, it is convenient to move to the interaction picture c ka (t) = e~ lH(,t c klJ (t)e ot . 
This transformation does not affect the correlation functions of our initial state (thermal equilibrium inside the source 
and vacuum outside), since it is invariant under the free evolution e~ lHot . We get 

(4 lffl (*)<*>»*(*)) = (ci iai (t)c k2<J2 (t)) 

P POC rOC 

^^\ 2 5 aia2 / d 3 kN{u k ) / dti [G(h)T]* klk e-^ / dt 2 [G(t 2 )T] k2k e^ . (4.6) 
J Jo Jo 

For small A and arbitrary t, one obtains (see Appendix |A"|) 

G kk ,(t) = S 3 (k - fc')e- I£fc * + 0(A 2 ) (4.7) 
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and in the weak-coupling regime we have 

(cUWcWt)) ^ ^ 1CT2 / d 3 fcAK) ^ in , Tk f in+ +Q(A 4 )- (4-8) 

j Efej - LUk + efe 2 — ^A; — iU^ 

All the other correlation functions are constructed from this two-point function, through the Wick theorem for an 
initial thermal state. 

It is instructive to write the correlation function (|4.8p in the configuration space: 

«(ri,t)^ a (r 2 ,t)) = S^p^ir^rz) \ 2 5 aia2 J d 3 fcJV(w fc )0j(ri)fo(r a ) + 0(A 4 ), (4.9) 

where 

fc ( r ) = / _£L Zk^ e lfc '- (4.10) 

is the Laplace transform of the free evolution of a wave packet 

r°° f d 3 k' i 

<Pk(r,s) = j o dtip k (r,t)e- st , <p k (r, t) = J -j== T k , k c l( - k (4.11) 

evaluated on the energy shell s = —xuJ k + + , i.e. ^fe(r) = (p k {r, — itJ k + + ). A particle with momentum k in the 
source is diffracted and propagates outside in the form of the wave packet ip k (r,t). The sum over k in formula (|4.9p 
yields the incoherent sum of such wave packets and a sort of "density matrix." 

V. BEAM PROFILE 

It is interesting to observe that the one-particle wave function (|4.10p can be expressed as a superposition of spherical 
wave originating from different points of the emitting region: 

&(r) = 0(ft,)/(-iV) J ^9(ro)v k °l (r), (5.1) 

with 

, . r H 3 k' 1 , „ife-r +ife|i — t | 

4°) (r) = / _jljL 1 lfe '.(r-ro) el fe-r = m ^ C (5 . 2 ) 

' J ^2^ iek ,-uj k -i0+ i\r-r \ V ' 

We intend to derive an expression which is valid far from the emitting region. Equation (|4.10|) reads 

Mr) = -^=- Tdpp 2 / d 2 p T{pp)k (5.3) 
V (2tt) 3 i Jo J Sp-Uk- 10+ 

For r — > oo, the phase p ■ r = cos 9 is stationary at 9 = and n, and the saddle-point approximation around these 
points yields 

k (r)~— ±=/ dpP 2 i(p#)fc / duuc-P™ 2 / 2 



'o 

/"OO 

J(-pf)fc lpr / duwe -p™V2 

£p - Wfc - l0+ 7 

1 r dpp T( ^ )fe e^, (5.4) 



which asymptotically behaves as 



mV2^9(k z )f(kr)g(Akr)- — , (5.5) 

lr 
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where (13. 5p is substituted for T p k and 

Akf = kf-k (5.6) 

represents the momentum transfer from k (before emission) to that directed towards position r with the same 
magnitude k (after emission). The Gaussian function 

g(Ak,) = JL_ e -^ w2 ^/ 2 (5.7) 

shows that particles with momentum fe in the source prefer to propagate in the same direction as k outside, but with 
some diffraction determined by the size of the window of the source. 



VI. CORRELATION FUNCTIONS 



We compute the spin-summed one- and two-particle distributions in the emitted beam, defined respectively by 

pP(r) = £ (ri(r,t)Mr,t)) = 2p?\r\r) (6.1) 

<r=T,J. 

and 

pf\r 1 ,r 2 ) = ^(n,^^,*)^^,*)^^!,*)) 

CT1,<7 2 =T,J. 

= ^(nlnJ^Cralra) - 2^ 1) (r 1 |r 2 )p| 1) (r 2 |r 1 ), (6.2) 

where p[^ (ri |r 2 ) was introduced in (|4.9[) . We are interested in the normalized two-particle distribution function with 
detector resolutions, 

n t= ~ \ P ( t 2 \r ll f 2 ) it(fi,f 2 ) ( . 

Ct(ri,r 2 ) = _ _ =1- _ _ (6.3) 



where 



Pt 1] (r) = J d 3 r-Rf=(r)p t (1) (r) = 2 J d 3 r Rf(r)p { t 1] (r\r), (6.4a) 
^(ri.fa) = J dV^fr,) / d 3 r 2 ^ 2 (r 2 )/3f ) (n,r 2 ), (6.4b) 



and 



^(n.fa) = 2 JdPnR^in) J d 3 r 2 ^ (r 2 )pf ) (n | r 2 )p^ 1) (r 2 |n ) (6.4c) 
arc defined in terms of the resolution function of the detector Rr(r), which is assumed to be Gaussian, 



o o d 

The quantity a characterizes the lateral size of the circular mouth of the detector and d the resolution in the longitu- 
dinal direction. The "interference term" 2t(fi,f 2 ) gives rise to a reduction in the two-particle distribution function, 
that is, antibunching. For bosons, the "— " sign in Eqs. (|6.2|) and (|6.3|) would be replaced by a "+" sign, and the 
coincidence count would be enhanced, exhibiting bunching. All the formulas below are easily switched to their bosonic 
counterparts by flipping the negative contribution of the interference term to a positive one. 
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A. Singlet and Triplet Contributions 

Before we compute the normalized two-particle distribution function (|6.3p . let us look at the structure of the 
two-particle distribution (|6.2p in the stationary beam: 

pp ) (r 1 ,r 2 )±^A 4 | d 3 fc x J d 3 k 2 N(u kl )N(cu k2 ) (4|£ fcl (n)| 2 |£ fca (r 2 )| 2 - 2^ 1 (r 1 )0 fcl (r 2 )^(r 2 )^ a (n)) 

= A 4 | d 3 fc! J d 3 fe 2 iV(a; fe JiV(a; fe j(3|<; i ) fc2 (r 1 ,r 2 )| 2 + |< > ) fe2 (r 1 ,r 2 )| 2 ), (6.6) 

where 

are the symmetrized/ antisymmetrized two-particle wave functions. Formula (|6.6[) for the two-particle distribution 
shows that 3/4 are contributed by the antisymmetric wave function while 1/4 by the symmetric one. This is because 
the thermal source is a complete mixture of the triplet and singlet spin states, the former being associated with an 
antisymmetric wave function in space, while the latter with the symmetric one, for the state of the fermions as a whole 
to be antisymmetric. A similar consideration applies to bosons, for which the symmetrized and antisymmetrized wave 
functions should be interchanged. 



^i(n) <£fei(f2) 



(6.7) 



B. Detector Resolution 



Let us now compute the normalized two-particle distribution function C s t{vi,V2) in (|6.3p in the stationary 
beam. The subscript "st" will henceforth denote quantities evaluated in the stationary limit, e.g. C s t(fi, r 2 ) = 
lim^oo Ct{f\, r 2 ). In the main part of this section, we shall employ an approximation which is nonsystematic but can 
nonetheless capture the essential features of the lateral effects of C s t ■ Its consistency and validity will be examined in 
Sec. EH] 

We need to evaluate the following component of the correlation functions in (|6.4p : by expanding r around the 
center of the detector, r = r + Sr, 

r ~ f + r ■ 5r 

r-f^ + lVf.Sr for Sr<f, (6.8) 

Akr ~ Akf. + *V*6r 
r 

with 7-V a projection operator which projects a vector v onto a perpendicular direction to f by V^v = v — f(f ■ v), 
we get 

d 3 r R,(r)^ ki (r)£ k2 (r) ~ ^^0(fc l2 )/( fcl #)0(fc 2z )/(fc 2 f ) C - Afc ^- w2A ^/ 2 e- A ^- w2A ^/ 2 c-'( fcl - fc2 ^ 



{2ir) 5 f 



x / (J 3 (57- p -8r-Bf^6r/2 p -(k 1 Ak 1 -+k 2 Ak 2 -)-W^V ii 6 r /f p -i(k 1 -k2)f-Sr 
^(27T) 3 dctP 2 



-i(fci— k 2 )r 

where 



A klk2 (f)Z klk2 (f)^c-^-^\ (6.9) 



A kl k 2 (r) 



(2tt) 5 ^/detp + (k'f + k'i)V 2 Vf,W 2 r^/f 2 } 



-Afc 1 --W 2 Afe lf /2 -Afe 2f -W 2 Afc 2 -/2 



x e (fciAfc 1 .+fc 2 Afc 2 .)-W 2 7^B 2 i 7^W 2 (fciAfc 1 .+^ 



Z fclfc2 (f) = ^(fci 2 )/(fc 1 f)^(/c 2z )/(fc 2 f)e-( fcl " fc2 ) 2f - B ^/ 2 , 



(6.10b) 
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and 

B^ 2 = V- 2 + (k 2 + k 2 )Vf.W 2 V^,/f 2 . (6.10c) 

Z klk2 (f) is responsible for the longitudinal effects and A klk2 (f) for the lateral effects. The one-particle distribution 
and the interference term in the two-particle distribution with detector resolutions are then given by 

pH\f) ~ 2A 2 1 J d 3 kN(uj k )A kk (f)Z kk (f) (6.11) 

and 

I s t(ri,r 2 )^2A 4 ^ f d 3 fcx / d 3 fc 2 N(u kl )N(u; k2 )A klk2 (n)^ lfc2 (f 2 )Z klk2 (rJZ^ (f^e" 1 ^-^-^. (6.12) 

r l r 2 J J 



C. Single-Particle Distribution 

Let us place our detectors on the longitudinal z axis, f = (0, 0, z). In this case, 



and therefore 



1 1 
I 1 

,0 0, 



Bi = 



l + a 2 w 2 {ki + k 2 2 )/z 2 






l + a A w 2 (k 2 1 +k 2 )/z 2 ° 





d 2 ) 



(6.13) 



A kl k 2 (r) 



(2w) 5 l + a 2 w 2 (k 2 + k 2 )/z 



_ e -w 2 kl_ L /2-wl(k 1 -ki I ) 2 /2 e -w 2 kl ± /2-wl{k 2 -k 2z ) 2 /2 



7 — J - J ^-^(k 1 k 1± +k 2 k 2± Y 



x e 



, (6.14a) 



Z klk2 (f) = e(k lz )f(k^)e(k 2z )f(k 2 f-)e-^- k ^ d I 2 . (6.14b) 
Now the single-particle distribution reads 

=a*_^ r^jfMi^f) r",.,,-^*-^..- (6 . 15) 



(2tt) 4 z 2 7 l + 2a 2 w 2 fc 2 /z 2 ./„ 

In order to estimate the integral over 9 by a Gaussian approximation, introduce a new integration variable by 

i + ww* **' + ro - t2(1 - cm8)2 ° G + ^ + ™ y ) ^ <6 ' 16) 

The above integral over 6* in (|6.15[) is reduced to the following form 

(p + q) f ' d9 p in( S ineco S e)-(p +(? ) s in 2 e (6.17a) 

Jo y p 2 cos 2 9 + q 2 sin 2 

P= l + 2 tt W/i 2 ' (07b) 

Notice that it is important to exponentiate all factors that change considerably in the integration region for the 
Gaussian approximation to be well posed. The exponent is expanded around its stationary point 

sin2e °^ 2(p + g) { 1+ P + 1- V± + (P + l) 2 ) (6-18) 
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and is approximated by a quadratic function of the form 
ln(sin 6 cos 6) - (p + q) sin 2 9 



\ + 2(vtqY ~) ~l( 1+P + q - v/l + b + 9) 2 ) - 2v/l + b + <7) 2 (e - 9 ) 2 . (6.19) 



2(p + g) 2 

The remaining slowly varying factor is estimated at 0q and we obtain 



' Vl + Cp + g) 2 -l c -i(i+ P + 9 -Vi+(p+9) 2 ) /■" /2 d9c -2 x /i+(P+9) 2 (e-e„) 2 

P 2 + <z 2 + (p - g)( v /l + (p + 9 ) 2 - 1) 



which, for large p> 1, is well approximated by 

Fk 1 Fk 



(6.20) 



-~7T\- for P »«■ ( 6 - 21 ) 



2yV + q/2 V e 2p V e 
Thus, the single-particle distribution (|6.15|) is evaluated as 



(2tt) 4 u; 2 z 2 y 7 ^2 cos 2@ + g 2 sin 2@ 



^ X \o \l~ I dfc7VK)/ 2 (fcf). (6.22) 

In the last line, it has been implicitly (and reasonably) assumed that the monochromator f(kf) extracts, in effect, 
only those momenta for which the inequality 

„..2j„2 

>max(l,iy z fc) (6.23) 



l + 2a 2 w 2 fc 2 /z 2 



holds. 

If the beam is well monochromatized around a given momentum fco and the distribution N(u)k) is a slowly varying 
function there, the one-particle distribution function (|6.22[) is further estimated for the Gaussian monochromator (|3.7[) 
as 



W - ^^^WL^W = V) 5 ^ 2 ^) 2 ^^" (6 ' 24) 



D. Two-Particle Correlation Function 



When the two detectors are placed on the z axis, i.e. n = r2 = (0,0, 1) = r, the interference term (|6 . 1 2[) reads 

r - ^ ^ , \ 4 2m 4 /■ 3 /■ 3 iV(^ 1 )iV(^ 2 )0(fc lz Mfc 2 ,)/ 2 (fc 1 f)/ 2 (fc 2 f) 
1st (n, r 2 ) ^ A Mn-9-9 / d fci / d fc 2 



where 



(2^) 10 z 2 z 2 7 7 [l + a 2 w 2 (fc 2 + fc 2 )/z 2 ][l + a 2 W 2 (fc 2 + fc 2 )/z 



2J 

1 -(fel-fc2) 2 d 2 n -l(fci-fc 2 )(2l-22) c -lU 2 fc? i -l«2( fe l~ fe l^) 2 c _to2 ' c 2_L~ ,U z( fc 2-fc2j) 2 



x e 



2 4 2 4 



2m fjs, frfi NMN^nkizMk^f^hr)/ 2 ^) 



= A* , in _,_ 9 d J fc x d^fc 



(2^) 10 z 2 z 2 i V 2 [l + a 2 «7 2 (fc 2 + fc 2 )/z 2 ][l + a 2 W 2 (fc 2 + fc 2 )/z 2 ] 

x c _ ^ fcl_fc2 ' 2d2 c~ 1 ' fcl_fe2 ^ ?1_?2 ^c _Plfe i i ~ 9l ' fcl_fcli ^ 2 e _P2fe ^^' j2 (' i:2_fc2 -) 2 e cfclJ -' fc2i 

(6.25) 

P < = U? - + a 2 l 2 Tfc 2 + fc 2 2 ) + If + JJ(k\ + kl) ) kl Qi = wl (l = 1 ' 2) (6 ' 26a) 
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and 



9 4 2 4 

aw aw 



z\ + a 2 w 2 {k 2 + k 2 ) + z 2 + a 2 w 2 (k 2 + k 2 ) J ^ (6 ' 26b) 

Integrations over the two azimuthal angles around the longitudinal axis yield a modified Bcsscl function of the first 
kind 

2tt ^ d0 12 c cfc - fc - cosfe = (2tt) 2 £ k\lkll = (27T) 2 I (ck 1± k 2± ). (6.27) 

The remaining integrations over k±± and k 2 ± can be performed just like before and we obtain [cf. (|6.15[) and (J6.21D ] 
2st(ri,r 2 ) 



(2tt) 8 z 2 z 2 J J 

y (c 2 /4)" n jVK,)/ 2 (fc t r) 



/ d0,- shift- e-P^ 8 ^-^ 1 - 



A 4 . 2 ^ 2 _ 2 /"dfcxfc? /^dfc 2 fc 2 2 e-( fel - fe2 ) 2d2 e-^- fc2 )^-^) 

(c 2 /4)" t r AT^J/^fc^) 



E l c TT 

,=o (rl!)2 ii 2 l + a 2 ^(fc? + ^)M 2 V d Pi J 2 Pl 



A 2(27r) 8 z 2 z 2 2 e J Q dfcl io dfc2 [1 + a 2 W 2 (fc 2 + fc 2 )/z 2 ][l + a 2 w 2 (k 2 + k 2 )/z 2 } PlP2 - c 2 /A 

= ^ -< r 4 4 2 2 - fV fdfa (6 . 28) 
2(2tt) s w 4 z 2 z 2 e J J 1 + a 2 w 2 (k 2 + k 2 ){l/2z 2 + l/2z 2 ) { ' 

For the well-monochromatized case, 

m 4 7T iV 2 K ) 



2st(rL,r 2 ) ~ A 4 



2(2tt) 8 u; 4 z 2 z 2 e 1 + a 2 w 2 k 2 {l/z\ + l/z 2 ) 

pOO p oc 

x / dfci / d/c 2 / 2 (fc 1 r)/ 2 (fc 2 f)e- (fcl ~ fc2)2d2 e" l(fcl - fc2)(Sl - S2) 



oo ./ — oo 

^ 1 exp(- T7 g^L |, (6 . 29 , 



2(2iry°w*z 2 1 z 2 (Sk ± ) i e 1 + a 2 w 2 k 2 (l/z 2 + l/z 2 ) ^/1 + A(6k z ) 2 d 2 *A l/(Sk z ) 2 + Ad 2 

and we end up with the analytical formula for the normalized two-particle distribution function 

— _ _ 1 1 1 / (zi - z 2 ) 2 

C st (r 1; r 2 ) - 1 - - i + a 2 w 2 k 2 {l/zf + l/z 2 ) ^ + A(Sk z ) 2 d 2 \ l/(Sk z ) 2 + Ad 2 > ' 



This is our central result. Its bosonic counterpart is readily obtained by just flipping the "— " sign in front of the 
second term. 

Before studying this formula numerically, it is worth analyzing its range of validity. 



E. Consistency of the Approximations 

One might wonder if the approximation and procedure we adopted when we performed the integration over r in 
(|6.9p are self-consistent, because, as some careful reader might have realized, we have partly kept second-order terms 
in Sr in the exponent of the integrand of (|6.9p . while only first-order corrections were considered in the expansions 
(|6.8p . Actually, we implicitly assumed that 5r = r — f is a small quantity, in order to keep only second-order terms 
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in 8r in the exponent, so that the integrals could be evaluated by Gaussian integrations. Stated differently, it is not 
clear whether we are allowed to expand the exponent of the integrand around f , since this is the stationary point of 
the exponent of Rr(r), but is not necessarily that of the integrand. In order for the approximation be consistent, we 
first have to find the true stationary or saddle point of the exponent of the integrand, r s , and then expand it around 
r s , keeping all second-order terms in r — r s . 

It is not difficult to derive the relation that the saddle point r s of the exponent of the integrand 



F(r) = -i(r - r) • 2?~ 2 (r - f) - i(fc x # - fe ) • (SK)- 2 (hf - fc ) - J(fc 2 r - fe ) ■ {5K)~ 2 {k 2 r - k Q ) 

- -Afc 1# • W 2 Afc 1# - ^Afc 2 f • W 2 Afc 2 r - i(fci - k 2 )r (6.31) 

has to satisfy. This reads 

= VF(r) = -V- 2 (r s - f) - h.p+ t (81C)- 2 far, - k ) - ^Vr s (5IC)- 2 (k 2 r s - k ) 
r B £r s Zr s 

- —Vf. W 2 Ak lf . - — Vf W 2 Ak 2f - i(fci - k 2 )f s . 



(6.32) 



The saddle point r s is therefore the solution of the equation 

r s =f- e(r B ), 

with 

e(r) = -V 2 Vfu(r) + i(ki - k 2 )V 2 r, 
r 



u{r) = (k 2 1 +k 2 )\W 2 + 1 -{5K.) 



r — 



ki + k 2 



(SIC)- 2 k -W 2 (k 1 k 1 +k 2 k 2 ). 



(6.33) 

(6.34a) 
(6.34b) 



This equation is iteratively solved to yield, after the first iteration, 



v — e(r) = z 



l-i(fe-fe)— # 



-(k\kx + k 2 k 2 ) 



(6.35) 



Notice that for this iterative solution to be a good approximation, its deviation from r must be a small quantity 
relative to z = |f|. Actually, this expression is still valid, after the second iteration, up to the second order in a/z and 
d/z. 

Since we obtain 



z 1 - \{k\ - k 2 )^r 



and r s 



-(fcifci + k 2 k 2 )_ 



(6.36) 



up to 0(a 2 /z 2 ) and 0(d 2 /z 2 ), the saddle-point value of the exponent, F(r s ), is easily estimated to be 
1,. . .o ,. . . 1 



F(r s ) = --(h - k 2 ) 2 d 2 - - k 2 )z - -w 2 {k 2 ± + fcfj 



2 A 



2z 



2 (fcifci + k 2 k 2 )\ 



(fci - feo) (k 2 ~k ) 2 1 2 2 2 . o r 7\ 



We then approximate the exponent by a quadratic form 



F(r) ~ F(r 8 ) + - ^(a* - a; si )fo - x sj ) q^qJ 



(6.38) 



to perform the Gaussian integration. The covariance matrix (the coefficient matrix of the quadratic term) reads 
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-a ( (z - z s ) + —(r- r s )± ■ b) - (r - r s )± ■ (m - — 1 (r - r 8 )_L, 



which yields, after integrations over r — r s 



a V det(.M - btf/a) 

Within the validity of the approximation adopted here, this factor is estimated to be 

a 2 w 2 {kl + ki) . i(fc!-fe 2 )a 2 



(6.39) 



(6.40) 



( , + + _ _ ) , ((U1) 



where 



A —1-1 1 1 A fc0(fcl+fc 2 ) \ ^| / fclfcl Z + fc 2 fc 2 A 

A fcl fc 2 1+ 2u;2((5fc±)2 2^2(^)2^ fc 2 +fe 2 J w 2^ k 2 + k 2 J' 

It is worth mentioning that we have obtained the same value as before [see, for example, (|6.14|) ] for the saddle-point 
value of the exponent (|6.37[) . within the validity of our approximation [up to 0(a 2 /z 2 ) and 0(d 2 /z 2 )]. The only 
correction to Ak 1 k 2 ( r ) in (|6. 14a|) is to replace its denominator by the quantity in the brackets in (|6.41l) 

a 2 w 2 (k 2 + k 2 ) a 2 w 2 (k 2 + k 2 ) ( , 2 \ 3d 2 - k 2 )a 2 

1 + _ 2 — » 1 + 4 ~ A fclfc2 + 2712— pT - — + - t. W ' 6 ' 43 

If we take, however, the same Gaussian approximation for the momentum integrations as in Scc. rvTCl by which the last 
term in the right hand side of (|6.42p is estimated to be ~ —w 2 /2w 4 (k 2 + k 2 ), and consider the well-monochromatized 
case, few corrections remain. Actually, it can be easily shown that the single-particle distribution (|6.24p has further to 
be divided by a factor 1 + (a 2 / z 2 )[k 2 / (5k±) 2 } and the denominator in the interference term (|6.29[) has to be corrected 
by 

1+aWfc °(w + w) ^ 1+aWfc Kw + w)( 1+ ^w)- (6 ' 44) 

It is remarkable and important to note that these corrections do not change the final result for the normalized two- 
particle distribution function (|6.30[) . illustrating the consistency and validity of the approximation and procedure we 
adopted there. 



F. Lateral Correlation at the Lowest Order for a Noncollinear Arrangement 

So far, we have investigated the lateral effects when the source and the two detectors are collinearly arranged. 
In order to bring antibunching to light, the second-order expansion of the exponent around the saddle point was 
necessary. It is interesting to note that, if the two detectors are placed off the longitudinal z axis, lateral ef- 
fects can be seen even in the lowest order with respect to 1/fi and l/f 2 . Indeed, when the detectors are placed 
on fi = (fi sin 8d cos $, f\ sin 8d sin fi cos 9d), r 2 = (— f 2 sin 9d cos $, — f 2 sin 9d sin <f>, f 2 cos 9d), the normalized 
two-particle distribution function (|6.3| is given by 



C (* *\ 1 ^i(Qd) / W 2 fc 2 (^) 2 sin 2 2e d (Sk z ) 2 (Sk 1 _) 2 (f 1 -f 2 ) 2 \ 

Cst(ri ' r2) = 1 ' 2 v /^2(9 d) ^3(9 d) CXP l~ 2D 1 (9 d )2 ?a (9 d ) ) ' (6 ' 45) 

provided the beam of particles is well-monochromatized and the two inequalities w^> w z , wk 3> 1 are satisfied. Here 
the auxiliary functions D\, D 2 , and D% arc defined as 

£>i(9 d ) = ( < 5fc 2 ) 2 sin 2 9 d + ((5fc_ L ) 2 cos 2 9d + 2w 2 ((5fc 2 ) 2 (Jfc_ L ) 2 (l-cos9d) 2 , (6.46a) 

D 2 (® d ) = D x (e d ) + 2w 2 (dk z ) 2 (6k ± ) 2 sin 2 9 d , (6.46b) 

D 3 (Gd) = D 2 (e d ) + 4a 2 (5k z ) 2 (Sk ± ) 2 sin 2 9 d + 4d 2 (5k z ) 2 (5k ± ) 2 cos 2 9 d . (6.46c) 
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As can clearly be seen, the two-particle distribution does depend on the lateral size w of the source. The derivation 
is discussed in Appendix iBl 

Now we consider two particular cases. In the first case where Gd = 0, one has -Di(O) = D 2 (0) = (Sk±) 2 and 
D 3 {0) = {Sk ± ) 2 [l + 4:(Sk z ) 2 d 2 }. Thus, (jtRS]) reduces to the lowest-order terms of ([00]) with respect to l/zi and l/z 2 : 

- ' - w+^r* ^'»/( w*) • <6 ' 47) 

where we have set z\ = f\ and z 2 = f 2 . In the second case, where r± = (x,y, z) and r 2 = (— x, —y, z), the two-particle 
distribution reads 

C st (^ 2 ) = 1 - T f ^W^+f)^ 2 ^ , ( 6 .48) 



2yjD 2 (x,y,z)D 3 (x,y,z) V D 1 (x,y,z)D 2 (x,y,z) 



where the functions D\, D 2 , and D$ are 

5i(x,y,«) = f^DiCGd) = (Sk z ) 2 (x 2 + y 2 ) + (Sk ± ) 2 z 2 + 2w 2 z (Sk z ) 2 (Sk ± ) 2 (f - z) 2 , (6.49a) 

D 2 (x,y, z) = f 2 D 2 {Q A ) = D x (x,y, 2) + 2 w 2 (<5fc z ) 2 (<5^) 2 (x 2 + y 2 ), (6.49b) 

5 3 (x,y,f)-f 2 i? 3 (e d ) = D 2 {x,y,z) + 4(Sk z ) 2 {5k ± ) 2 [a 2 (x 2 + f) + d 2 z% (6.49c) 

and f = \/ x 2 + y 2 + z 2 . The geometrical configuration of noncollinear detectors is not the usual one and will not be 
discussed any further. It might be important for electron antibunching from superconducting emitters [23| . 



VII. LONGITUDINAL, LATERAL, AND TEMPERATURE EFFECTS ON ANTIBUNCHING 

The normalized two-particle distribution function C B t(fi,f 2 ), evaluated on the basis of the expressions (|6.22|) and 
(|6.28|) . is shown in Figs. [2HU The two-particle distribution (the number of coincidence counts) is always suppressed 
when the two detectors are close together. The dips in the figures represent antibunching. 

It is clear from the expression (|6.2j) for the two-particle distribution that the number of coincidences is reduced 
to one half of that naively expected on the basis of the counts by the single detectors, when the two (ideal) point 
detectors are at the same point. This is understood by the expression (|6 . 6|) : the triplet spin states accompany the 
antisymmetric wave function in space, yielding antibunching, while the singlet spin state accompanies the symmetric 
wave function in space, yielding bunching. The interplay of these contributions (three fourths from the antibunching 
and one fourth from bunching) results in the minimum value 0.5 (= 1 — 3/4 + 1/4) of the normalized two-particle 
distribution function. The width of the dip, on the other hand, is governed by the width of the spectrum of the 
emitted particles, as is clear from the analytical formula (|6.30p or from the Fourier-integral representation of the 
interference term in (|6.12[) . 

Let us next look at the effects of the detector resolutions. Not only the longitudinal resolution of the detectors d 
(Fig. [2]) but also the lateral size of the detectors a affects the visibility of the antibunching (Fig. [3]), although we are 
looking at the coincidences between the two detectors located along the longitudinal axis. The width of the dip is 
broadened by the longitudinal resolution of the detectors d, while it is not by the lateral size a. 

The analytical formula (|6.30p clearly reveals how the resolutions of the detectors affect the coincidence counts. It 
also shows that the detectors can be regarded as point detectors when 

a « J - , rf «^TT ; 
V2wk 2dk z 

and these quantities characterize the lateral and longitudinal coherence lengths, respectively. Clearly, the above 
conditions agree with those derived in classical textbooks 0| . 

The temperature of the source affects antibunching in a way that deserves a few words of explanation. The visibility 
(namely, the depth of the dip) is temperature independent, as a consequence of the antisymmetry of the fermionic 
state a la Pauli's principle: indeed, antisymmetry is exact, both for pure and mixed states, and is preserved even at 
very high temperatures, where the fermionic state is totally mixed. See Fig. Sfa). This is clear in our formulas and 
figures: the prefactor —1/2 does not depend on any details of the source and the experimental setup. See e.g. Eq. 
(|6.6p : in the second line, irrespectively of the temperature distribution N(u>k), when |n — r 2 \ — > oo the interference 
term (second term in brackets) goes to zero, while at T\ = r 2 it equals half of the background term (first term in 
brackets). On the other hand, the width of the dip depends on temperature and even strongly by the location of the 
monochromator window with respect to the Fermi level, as can be clearly seen in Fig. (Hb). As a consequence, the 
dip becomes narrower as temperature is increased and the effects of antibunching becomes more difficult to detect. 
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FIG. 2: (Color online) Normalized two-particle distribution function C s t(0, 0, fi; 0, 0, 160k;) vs the longitudinal detector coordi- 
nate z\ . Cst is evaluated in the stationary state on the basis of the expressions (|6.22[) and (|6.28[) , with the Gaussian detectors 
located at f\ — (0,0, zi) and f"2 = (0,0, 160w). We focus here on the effect of the longitudinal resolution of the detector, d. 
The parameters are fco = 20 w~ , 5k z = 0.5 w~ , a = 0, and from bottom to top (in the dip) d/w = 0, 1, 2, 3, 4, 5 (in units h — 1 
and m = 1). We set /3 = 5 rnw 2 /h 2 (a low temperature) and the Fermi level fi — (ko + 8k z ) 2 /2m ~ 210 ti 2 /mw 2 (just above the 
momentum window). Note that in this case p ~ 400 in (|6.17p and the condition (|6.23p imposes w z <C 20 w. The values based 
on the numerical integrations of (|6.15[) and (|6.25|l without the Gaussian approximation are also shown by dots for d = and 
• r )w;. These were checked to be independent of w- for small w z . 



It is useful to summarize the meaning of our analysis. Equations (|6.ip and (|6.2j) yield the one- and two-particle 
distributions in the beam, that are expressed in terms of the fermionic operators. Equations ()4.9j) and (|6.6p then 
express these quantities in terms of the "wave functions" (f>k(f) and of the temperature dependent function N(loi;). If 
the Fermi distribution is plugged in, all formulas apply to fermions: otherwise the analysis in Sees. ITVl lVIII is general 
(modulo some sign changes) and can be applied to bosons as well. 

It is interesting to apply our final result (|6.30[) to some interesting experimental situations. It is necessary to stress 
that our analysis is strictly valid only for experiments such that the beam of emitted particles travels in vacuum. If this 
situation closely resembles the experimental one, then our equations apply; otherwise, additional care is required in 
order to explain the experimental data. In some experiments, like those in which correlation in the current intensities 
are observed fiol . fllT | , our formulas cannot be applied and a different analysis is required (23j . 

Let us start from an analysis of the electron experiment [l2l ]. One infers the values a ~ 2mm, w ~ 18 nm, 
fco — 10 1:1 ni _1 , z ~ 10cm. By plugging these values in Eq. (|6.30p . one sees that the first of the two factors 
multiplying the exponential is of order 10 -4 . Moreover, the coherence time is Ai co h — 32 fs, while the response 
time of the detectors Aidct — 26 ps, which yield a value ~ 1/300 for the second factor in front of the exponential. The 
global factor multiplying the exponential is therefore of order 10 -6 , which makes the observation of the phenomenon 
quite complicated. Indeed, the authors had to apply a lateral magnification technique (nominally of order > 10 4 ) in 
order to observe antibunching. Notice also that in our formulas the Coulomb repulsion is neglected. This is a delicate 
issue that would require additional investigation. 

Let us now look at the neutron experiment [l7| . The relevant values are a ~ 1 cm, w ~ 1 cm (a mosaic crystal 
was used in order to reflect the beam into the apparatus), ko ~ 10 10 m _1 , z ~ 10m. The beam coming out of an 
oven travels in waveguides for about 100 m, is then monochromatized through back scattering by a prefect crystal 
and illuminates the whole mosaic crystal on a region of order few cm 2 , the back reflection being coherent only on 
regions of order ~ fim 2 that are uniformly distributed in the whole volume. By plugging the numerical values of the 
parameters in Eq. (|6.30[) . the first factor is of order 10 -10 . Moreover, by comparing the coherence time of the neutron 
wave packet A£ co h ~ 20 ns with the response times of the detectors Aidct — 0.1 /is (two different types of detectors 
were used, with response times that differ by a factor 10-20), one obtains a second factor of order > 10 _1 (or smaller 
by a factor 10 for the other type of detectors), yielding a very small antibunching dip. It is interesting to observe 
that if we take w ~ 1 /mi (the size of a monocrystal in the mosaic), by Eq. (|6.30|) the first factor is of order 10 -2 , the 
second factor remains identical and one obtains an antibunching dip of a few percent, which can be brought to light 
by deconvolution and is in agreement with the experimental data. An exhaustive analysis of the physical effects of 
the mosaic crystal used in back reflection is involved and will be presented elsewhere. 



Another interesting experiment is that performed with X-rays 12 1| . An important part of this experiment is devoted 
to the analysis of the lateral coherence features of the beam. The setup involves the values k ~ 7 • 10 10 m _1 , z ~ 70 m, 



VIII. APPLICATION TO EXPERIMENTS 
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FIG. 3: (Color online) (a) Same as in Fig. [2] but with d — and from bottom to top a/w = 0, 1,2,3,4,5. We study here 
the effect of the lateral size of the detector, a. Note that p ~ 225 for a = 5w in (|6,17[1 and the condition (|6.23[) imposes 
iiij <C 11 to. The values based on the numerical integrations of (|6.15[) and (|6.25|) without the Gaussian approximation are also 
shown by dots for a = and 5w and were checked to be independent of w z for small w z . (b) Depth of the antibunching dip, 
1 — C s t{0, 0, 160w; 0, 0, 160w), as a function of a. All the parameters are the same as in (a). 



0.8 - 



o 0.6 



0.4 



'O 0.2 - 




0.0 



154 156 158 160 162 164 166 

Zi/w 



1.2 
1.0 

0.8 
0.6 
0.4 
0.2 
0.0 









(b) 











4 6 8 
wk, wk z 



10 



FIG. 4: (Color online) (a) Same as in Fig. [2] but with a = 0, d = 0, and from bottom to top f3/(mw 2 /h 2 ) = 5,0.2,0.05. We 
analyze here the effect of the temperature, /9" 1 . The values based on the numerical integrations of (|6.15p and (|6.25[) without 
the Gaussian approximation are also shown for (3 /{mw 2 /h 2 ) = 5 and 0.05 by dots, and were checked to be independent of 
w z for small w z . (b) Particle emission (from the oven into the beam) is proportional to the overlap of the form factor of the 
monochromator (bell-shaped curve) and the Fermi distribution function in the oven (step- like function). The effective width 
of the overlap is inversely proportional to the width of the (antibunching) dip in (a). 



while the detector size was changed in the range a ~ 10-400 /iin. The authors observed a reduction of the bunching 
peak very similar to that described in our Fig. [3Jb), as a function of the size of the detector mouth a, and estimated 
the lateral coherence length of the beam from the width of the plot, obtaining a source size of order w ~ 20 ^m. The 
longitudinal coherence length is known to be A£ co h — 2 mm, while the length of the photon bunch is A^dct — 5 mm: 
by plugging these values in our Eq. (|6.30[) . the first factor is ~ 0.9 for a ~ 10 fxm, while the second one is ~ 0.4, so 
that (1/2) ■ 0.4 • 0.9 ~ 0.18, which explain well the observed data (a positive bump ~ 0.25.) 

Finally, it is interesting to look at some experiments done with thermal optical photon, in order to study both their 
correlations and imaging. Experiments of this kind were widely debated during the last few years @, H, 0, We 
shall focus on the experiment (25| , in which the source is a He-Ne laser (in other similar experiments a green-doubled 
Nd:YAG laser was used). One estimates a ~ 30 /im, w ~ 1mm, ko ~ 10 7 m~ 1 , z ~ 150mm. The pseudothermal 
source is obtained by randomizing the phase of the photon field by means of a rotating ground-glass disk (so that the 
expression random source would probably be more appropriate). The first factor in Eq. (|6.30|) is therefore ~ 0.1 while 
the second one is essentially 1 (Atdot — 500 ps, yielding A^dct — 0.15 m, which is to be compared to the much larger 
coherence length of the laser). This yields an overall dip of good visibility. Unfortunately a quantitative comparison 
with the experimental data is difficult because the authors, being interested in the change of the width of the second 
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order correlation function with the source size, only plotted a (re) normalized correlation function. 

This brief summary of experimental data shows that our analysis and final formulas agree well with most experiments 
performed so far, with different particles. In some cases a comparison is more complicated and/or requires additional 
information. One phenomenon that we find of interest, but still lacks experimental confirmation, is our prediction 
that for fermionic systems we expect that the visibility of second-order interference effects should show no dependence 
on temperature, as explained in Fig. [4] We emphasized at the end of Sec. I VIII that this is due to the exactness of 
Pauli's principle (yielding perfect antisymmctrization) even for mixed states. 

IX. CONCLUSIONS, COMMENTS, AND PERSPECTIVES 

We analyzed antibunching in a beam of non-interacting fermions and investigated the behavior of the visibility as 
a function of the size of the source and the detectors, as well as the temperature of the source. These parameters are 
critical and play a prominent role in experimental applications. Our analysis makes use of Gaussian functions both 
for the emitting region of the source and the detector, and is adapted to an approximately cylindrically symmetric 
situation, with circular detector placed close to the longitudinal axis. This is the relevant situation in most experiments, 
in particular with neutrons and electrons (where however additional Coulomb effects, as well as more specific emission 
features of the source need to be scrutinized). It is also worth noticing that in the observation of pion correlations 
fl6l ] the experimental data have been exploited to determine the dimension and the expansion dynamics of the pion 
"source" (fireball produced in central Pb-Pb collisions, which is expected to be a droplet of quark gluon plasma at 
the freeze-out point). Clearly, our approach can prove to be quite useful in such a situation as well, when source size 1 
and temperature are not known. 

Let us look at possible applications and future perspectives. First of all, we emphasize that although our analysis 
was performed for fermions, all formulas can be easily translated to the case of bosons, enabling one to scrutinize 
interesting novel experiments in quantum imaging and lithography. Some recent applications make use of chaotic or 
pseudo-thermal light sources, to which our formalism immediately applies 0, Br- 
other possible applications are in solid state physics, where, as a consequences of the symmetrization procedure of 
the many-body wave function, entanglement should be present in bulk matter, raising delicate problems in relation to 
its detection and extraction p6l . |27l |28| . Since one should get entangled neutron pairs within the coherence volume of 
the wave packet (antibunching being observed within the same coherence volume), these pairs could be used as very 
efficient "probes" for entanglement in solids in future experiments. This is clearly relevant for quantum information 
processing and for tests of the Bell inequalities. 

Finally, we mention some interesting speculations in neutrino physics and the structure of the universe, where, 
as a consequence of entanglement formation, the hypothesis of a collisionless fluid of classical point particles can 
be critically re-examined, yielding a "quantum overpressure," with significant consequences during the non-linear 
structure formation epoch at low redshifts [29| . It is remarkable that independent ideas can bear consequences in very 
diverse fields. 
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APPENDIX A: G(t) IN THE WEAK-COUPLING REGIME 

Let us prove (|4.7[) . We write the matrix (|4.5b|) as 

G- 1 (s) = D(s) + X 2 K(s), D kk ,(s) = (s + ie k )6 3 (k-k'). (Al) 

Let us show that 

[l + X^K^D-HsT 1 (A2) 
docs not have poles in the weak-coupling regime, i.e. 

Det[l + X 2 K(s)D- 1 (s)} ^ 0, (A3) 
for a sufficiently small A. Indeed, the determinant is evaluated as 

Det[l + \ 2 &(s)£>- l (s)] = c ^o g [i+x 2 k( s) D-\ s )]^ (A4) 

which reads, in the weak-coupling regime, 

= l + A 2 Tr[A'( S )i)- 1 ( S )]+0(A 4 ) 

= 1 + A 2 / d 3 fc^^ + 0(A 4 ) 
J s + ie k 

= 1 + A 2 /d 3 fc /d 3 fc' 7 r + 0(A 4 ). (A5) 

J J (s + ie k )(s + iuj k >) 

Let us look for a zero of the determinant on the first Riemannian sheet, by putting s = —iw — 7/2: 

DC[1 + - 7 /2)] . 1 - A* J A J *V fo _, + , 7 ^_„ + 17/2) + 0(X'). (A6) 

Assume now that the square of the emission matrix in the energy representation, 

r(w, u/) = (2tt) 2 / d 3 fc / d 3 fe' \T kk , \ 2 5{e k - u)S(u k > - to'), (A7) 



docs not contain any "nonlocal" part 5(u> — u>') and has the properties 

r(w,u/)^0 for w,o/->-0,oc (A8) 

[assuming some good continuity property for T(ui,ui')]. For such a "reasonable" emission matrix T kk i, the integral in 
the second term of (|A7j) 

dbJi [°° du 2 r(wi,w 2 ) ^ AQ ^ 



o 27r J 2tt (uji — cj + 17/2) (<J2 — u + 17/2) 

is convergent for any u> and 7, and the determinant can always be made non-zero by choosing a sufficiently small A. 
The pole of 

G(s) = D-\s)[l + X 2 K(s)D- 1 (s)]-\ (A10) 

therefore, comes from the first factor D^,(s), i.e. s po i c = — i£fc, and we are allowed to expand the nonsingular second 
factor as a power series of A, yielding 

G kk ,(t) =S 3 (k-k')e- ls fct + 0(A 2 ). (All) 
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APPENDIX B: TWO-PARTICLE DISTRIBUTION IN A NONCOLLINEAR ARRANGEMENT 

Here we briefly sketch the derivation of (|6.45|) . Remembering 

fi = (fi sin Qd cos <£>, fi sin 8d sin i>, fi cos 0d), (Bla) 
f-2 = (— f2 sinOd cos$, — f 2 sin9d sin<I>, f 2 cos 6d), (Bib) 

the functions Z klk2 (f) and ^4fc 1 fc 2 (T 5 ) read 

Z klk2 (f) = 0(fc l2 )/(fc 1 f)0(fc 2z )/(fc 2 f)c-( fcl - fc2 ) 2f ' I,2f / 2 , (B2a) 

^ felfe2 (f) = ^L c -(Afe 1 ,-W 2 Afc lf + A fc2f .W 2 A fc2 ,)/2^ (B2b) 
(27TJ 5 



By substituting (|B2[) into (|6.1ip and (|6. 12[) and introducing spherical coordinates, we get 

Pst(r) = A 2 -^-j / dkk 2 N{u k )f 2 {kf) MsmO d<£e- Afe * w2Afe *, (B3) 
l/7Tj r Jo Jo Jo 

J s t(n,f 2 ) = A 4 _ 2 _ 2 / dfc!fc 2 / dk 2 klN(Lj kl )N(Lj k2 )f(k 1 h)f(k 2 h)f(kih)f(k2f-2) 
{^n) r 1 r 2 j Jo 

where 

J(fe;ri,r=a)= dSmnO (** d^-CA^-w^+A^.wM^)^ (B5) 
Jo Jo 



Note that, when ri and r 2 satisfy ()B1[) . one has 

. . 1 / fc 2 sin 2 6 d (fccos9 d - A: ) 2 \ 
/(fcn) = /(fcr 2 ) = [(2 , )3 (^ )2( ^ ) 4 ]1/4 exp( v - 1 ^ IF - 4((5fcz)2 ) - /(fc.Od), (B6) 



and fi • £> 2 fi = f 2 • £> 2 r 2 = a 2 sin 2 9 d + d 2 cos 2 6 



d- 



1. First-Order Correlation 

Let us first consider the case where r = r%. Since 

Afe # • W 2 Afc # = w 2 fc 2 (sin 2 e d + sin 2 6») + w 2 fc 2 (cos6 d - cos(9) 2 - 2w 2 k 2 sm6 d sin(9cos($ - 0), (B7) 

one has 

rir/2 p2n fir/2 r2ir 

dOsmO d(j>e~ Ak r w2Ak r =e - psin2ed / d6»sin6»e- psin2e - <?(cosed - cose)2 / d0e 2psinedSin9cos (' 1 '-^, 
Jo Jo Jo 

(B8) 

where p = w 2 k 2 and q = w 2 k 2 . The integral over (j> yields a modified Bessel function of the first kind: 

/ 2T d0e 2 PS ine dS in e co S (*-0) = 27r / (2p sin d sin 0) = 2tt V (P sh ^^°) 2n \ (B9) 
Jo „^o ( m! ) 



Therefore, with the aid of the trick used in (|6.28p , we have 

tt/2 p2n 

d0sin0 / d<^»e- Afe * w2Afe * 

(p 2 sin 2 e d ) m r /2 



o Jo 



oo / 2 • 2f\ \m pn/2 
27re -psin 2 e d \P Sm Vd) / d6)gin 2m+l 6)e -psin 2 (?-< ? (co S e d -co S e) 2 

(^O 2 Jo 
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oo 



[m'.y \ op I I- 



m=0 

Hereafter, the integral over 9 is evaluated when p 3> q and p> 1. By changing the integration variable, we have 

tt/2 -P( 1+ 'T-T d ) / j ^z^+^^fi \ 

d6»sin6»e- psin2e - 9 ( cosed - cos ^ 2 = / dxe* + / * " dxe x * . (Bll) 



Vp-i 



The first term is of order of e p / y/p. Integrating by parts, the second term becomes 

*- 2 M „„„„«. » / \ 

l + 



-PU+-S35- 4 ) „./p=?+^§^ 2 e - 9 (co S e d -l) 

da; e^ 



\/p~~q J i 2(p - g + gcos9 d ) 



f d^, (B12, 



3ee V J c \ / fV¥=g+^rf 3e 



4 Vp-g VP - g 

where the second and third terms are, respectively, 0(l/p 3 ^ 2 ) and 0(e~ p /y/p), and the last term can be shown to be 
0(p" 5 / 4 ). Thus, one has 

p*/ 2 , , „-?(cose d -i) 2 

/ dflsin 9 c- psin e-9(™ s e d - cos e) = + (p-s/4) 

7o 2(p-g + g cose d ) 

-«j(cos e d -i) 2 

= Y p + °( p " 5/4 )- ( B13 ) 

By substituting (|B13[) into (|B10[) and retaining the leading order terms in 1/p, we have 



/ d0sin0 / d0e- A ^ w2Afc * = 2 TC - psin2ed V ■ 

JO J m =0 



(m!) 2 \ dp 



sin 2 e dp -?(cose d -i) 2 (p 2 sin 2 e d ) m m! 

/ 4 (m') 2 



7re H a e 



m=0 

^_ — psin 2 © d — q(l — cos d ) 2 psin 2 d 
P 



Jl_ -to 2 fc 2 (i-co S e d ) 2 



l k 2 



(B14) 



Then, the one-particle correlation is 



^t(r) = A 2 2 _ 2 / dkN(co k )f(k,Q d )e-^ k ( 1 ~ cose ^ . (B15) 

When f = f2i we obtain the same result. This is the denominator of Eqs. (|6.3j) and (j6.45j) . 

2. Second-Order Correlation 

Because (|B1[) leads to 

i(Afc fi •W 2 A% 1 +Afc #2 •W 2 A% 2 )=p(sin 2 ed + sin 2 0) + g (cose d -cos6l) 2 , (B16) 
the auxiliary function J(k;fi,f2) is evaluated as 

f-n/2 



J(k; fi, T2) = 27r / d^sin^e 
Jo 



7T _,„2j,2 =; „2 

w 2 fc 2 



p(sin G d +sin 0) — g(cos d — cos f 
uj^fc* sin i! e d „-Mi 2 fe 2 (l-cose d ) 2 



(B17) 
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where we have used (|B13|) in the second equality. Thus, in terms of f(kj, O d ) (j = 1, 2), 



X st (r 1 ,r 2 ) = A 4 



OO pOO 



2(2n)Wf 2 f 2 j q j Q 
This is the numerator of Eqs. (|6.3| and (|6.45l) 



dfc! / dfc 2 iV( Wfel )^(u; fe2 )/ 2 (fc 1 ,e d )/ 2 (/c 2 ,e d )e- I ( fcl - fc2 )^-^) 



x e -(fc 1 -fe 2 ) 2 (a 2 sin 2 e d +d 2 cos 2 e d ) e - (fc 2 +fc 2 ) [w 2 sin 2 e d +to 2 (1-cos e d ) 2 ] (B18) 



3. Well-Monochromatized Case 

If the beam of particles is wcll-monochromatized and the distribution N(ujk) is a slowly varying function there, we 
have 



Art(r) = A* 2 _ 2 iVK n ) / dfc/ 2 (A; ! e d )e-^ fc2 ( 1 - cose ^ 2 



{2w) i w 2 r' J _ 00 

9 

= A 2 - — - — . , . iV(u; fco )exp 



(2 7 r)5 M ,2r 2 (,5fc J .) v AD i (e d ) 



fa 



2(5k z y 



1 



(a-_L) 2 cos 2 8 d 

r>i(e d ) 



(B19) 



where the Gaussian fc-integration has been carried out with the aid of 

w 2 fc 2 (l - cos9 d ) 2 



k 2 sin 2 9 d k 2 cos 2 9 d — 2k k cos d + k 2 



2(Sk ± y 



On the other hand, we have 



2(5k z y 



k cos 6 d k 2 

- rC 



2(8k z ) 2 (6k ± ) 2 (6k z ) 2 2(Sk z y 



(B20) 



Z s t(ri,r 2 ) = A 4 



2(2 7 r) 8 w 4 f 2 f 2 / 



OO pOO 



OO «/ — oo 



N 2 (u; ko ) / d*i / dfc 2 / 2 (fc 1 ,e d )/ 2 (fc 2 ,9 d )e- I ( fcl - fc2 )^- F2 ' 



x e 



(fcl-fc 2 ) (a 2 sin 2 e d +d 2 cos 2 e d ) -(k(+k%)[w 2 sin 2 e d +u) 2 (1-cos 6 d ) 2 ] 



(B21) 



In terms of K = (k\ + fa 2 )/2 and fa = fci — fc 2 , one has 



2 /fc 2 sin 2 8 d fc 2 cos 2 e d -2fc fcjcose d + fc 2 \ „ iw _ 

g l^^piF + W J " 1(1 " 2)(fl " f2) 

- (fci - fc 2 ) 2 (a 2 sin 2 9 d + d 2 cos 2 9 d ) - (fc 2 + k 2 )[w 2 sin 2 9 d + w 2 (l - cos 6 d ) 2 ] 



£>2(©d) r ,2 2fc COs9 d 



A' 



£> 3 (e d ) 



(<5fc z ) 2 (<5fciJ 2 (<5fc z ) 2 (<5fc 2 ) 2 A(Sk z ) 2 (5k ± ) 2 



: k 2 - ik(n - fa), (B22) 



and, thus, 
2st(n,r 2 ) = A 4 



Tft 



2(2Try°w*f 2 f 2 (6k ± ) 2 ^D 2 (e d )D 3 (e d ) 



x exp 



/.• 2 



1 - 



(<5fc J .) 2 cos 2 e d \ (5fc z ) 2 (5fc J .) 2 (f 1 -f 2 ) 2 



(<5fc,) 2 D 2 (e d ) ) D 3 (e d 

By plugging (|B15[) and (|B18|) into (|6.3[) . we obtain the normalized two- particle distribution function (|6.45[) 



(B23) 
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